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We study intrinsic simulations between cellular automata and introduce a new necessary condition for a CA to simu- 
late another one. Although expressed for general CA, this condition is targeted towards surjective CA and especially 
linear ones. Following the approach introduced by the first author in an earlier paper, we develop proof techniques to 
tell whether some linear CA can simulate another linear CA. Besides rigorous proofs, the necessary condition for the 
simulation to occur can be heuristically checked via simple observations of typical space-time diagrams generated 
from finite configurations. As an illustration, we give an example of linear reversible CA which cannot simulate the 
identity and which is 'time-asymmetric', i.e. which can neither simulate its own inverse, nor the mirror of its own 
inverse. 
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1 Introduction and definitions 

Cellular automata (CA) are well-known for the variety of behaviors they can exhibit. A lot of classifica- 
tion schemes where proposed in the literature, trying to make this variety of behaviors more intelligible 
MWol84 , Gil87 , Kur97 1. Such classifications usually consist in a (finite) list of distinctive properties giving 
raise to a partition of the class of all CA. Another approach consists in defining a simulation relation be- 
tween CA, and studying the ordered structure induced by the simulation. We follow this latter approach, 
and more precisely the simulation relation defined in I D MOT 1 1 allDMOT 1 lb I giving rise to the notion 
of intrinsic universality [O1108]. The intuition behind this simulation relation is simple: a CA is simulated 
by another if some rescaling of the first is a sub-automaton of a rescaling of the second. 

More formally, we restrict ourselves to dimension 1 and the definition is as follows. A CA F is a 
sub-automaton of a CA G, denoted F C G, if there is an injective map ip from A to B (state sets of F 
and G respectively) such that Ip o F = G olp, where Tp : A z — > B 1 denotes the uniform extension of ip 
to configurations. We sometimes write F C y G to make tp explicit. This definition is standard but yields 
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a very limited notion of simulation: a given CA can only admit a finite set of (non-isomorphic) CA as 
sub-automata. Therefore, following works of J. Mazoyer, I. Rapaport and N. Ollinger |MR98 O1102 
DMO TTTal IDMOT1 lbl . we will add rescaling operations to the notion of simulation. The ingredients 
of rescaling operations are simple: packing cells into blocks, iterating the rule and composing with a 
translation (formally, we use shift CA a z , z € Z, whose global rule is given by a z (c) x — c x _ z for all 
x E Z). Given any state set Q and any m > 1, we define the bijective packing map b m : Q z —> (Q m ) Z 
t\v: 

Vz £ Z : (& m (c)) (z) = (c(mz), . . . , c(mz + m — 1)) 

for all c £ Q z . The rescaling F< m >t> z > of F by parameters m (packing), t > 1 (iterating) and z G 1 
(shifting) is the CA of state set Q m and global rule: 

b m oa z o F 1 o &" 1 . 

With these definitions, we say that F simulates G, denoted G =4 F, if there are rescaling parameters mi, 
?7i2, h, t 2 , z 1 and z 2 such that G <mi ' tl ' Zl> C F <m2 < t2 ' Z2> . 

Determining whether some given CA simulates another given CA is hard (undecidable in general 
IDMO TTTbl section 4.3]). For instance, looking at typical space time diagrams of two CA gives no clue 
on whether one simulates another, because the simulation can occur on a set of configurations of measure 
0. Despite the general undecidability of the simulation relation, one can still hope to better understand 
its restriction to some specific classes of CA. For instance, the simulation relation is fully understood 
on products of shifts ID MOT1 lbl theorem 3.4] thanks to a 'characteristic sequence' which is essentially 
the sequence of ratio of translation vectors. Hence, if F = Co x a\ x 03, one can prove that F cannot 
simulate F^ 1 = ctq x er_i x cj_3 because they do not have the same characteristic sequence. 

In this paper, we introduce a general necessary condition for a simulation between two CA to be possi- 
ble. It focuses on surjective CA, but we will essentially use it on linear reversible CA. This condition is 
expressed as a characteristic set \ °f points of the real half-plane which is decreasing w.r.t. =<! (Theorem[T| 
below): 

F^G^ X (G)C X (F). 

A striking property of \ is that it can be somewhat visualized on typical space time diagrams of linear 
CA. Moreover, the set \ is closely related to so-called 'Green functions' of linear CA for which sys- 
tematic analysis techniques have been developed in HGNWIOL Hence, formal proofs of impossibility of 
simulation between two linear CA can be derived from heuristic observations of space-time diagrams in a 
quasi-automatic way. 

The set of reversible CA is somewhat structured with respect to =3! since it possesses a maximal element 
(i.e., a reversible universal CA ID MOT1 lbl theorem 4.5]) and verifies the following ID MOTlTbl theorem 
4.4]: 

F^G^F- 1 ^ G- 1 

Therefore, a reversible CA is either ^-equivalent to its inverse, or ^-incomparable to it. The most complex 
reversible CA, reversible universal CA, are all ^-equivalent to their own inverse. Coming back to the 
example F above (product of shifts), we have that F and F^ 1 are ^-incomparable. Following MAN 101 . 
let us associate to every reversible CA F its dual F — A/oF^'o M, where M is the mirror transformation 
on configurations (M(c) z = c_ 2 ). Any product of shifts is self-dual, and generally speaking it seems to 
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be hard to come up with CA that do not simulate their dual, while non-time-symmetric CA in the sense of 
|MG10| come in profusion. 

An interesting question in this context is how different a reversible CA can be from its dual. As an 
illustration of the necessary condition for simulation between CA that is given by Theorem [T] we study 
in section [3] some reversible linear CA. The first one simulates its inverse, its mirror, its dual, but not the 
identity; the second one simulates neither the identity nor its inverse or its mirror image or its dual. 

2 Simulation and geometry 

The basic ingredient in this section is the collection of functions telling how a change of value of the 
center cell in the initial configuration will affect some other cell's value at some step in the future. Such 
functions are often studied for linear cellular automata (see section [3) and are sometimes called 'Green 
functions' in this context |Moo98|. 

Let F be any CA and fix some x £ Z and some y £ N. For any configuration c £ Q z and any q £ Q, 
we denote by (f> c (q) the following configuration: 



Mi)?. = 



q if z = 0, 
c(z) else. 



We then denote by F| c : Q ->• Q the map q i-» (Fy{cj> c {q))) x . 

For instance, if F is simply the identity, F$ c is the identity when x = 0, otherwise it is the constant 
function q H> c(x). For a less trivial example, consider the cas where F = (J) is the sum with neigh- 
borhood {0, —1} over Q = Z/2Z, i.e. @(c) x = c(x) + c(x — 1). Starting from a single nonzero cell, 
iterations of this automaton generate Pascal's triangle modulo 2. For x £ N, let x = Yl n & b x (n)2 n , with 
b x (n) £ {0, 1}, be its binary representation, and B x — {n £ N\b x (n) = 1}. Then 



V ( \ _ / ® V { c )x + q if x > and B x C B y 
*,c {<1) ~ { @ v (c) x else 



We are interested in positions in space-time where the influence of the center cell is concentrated, 
whatever the initial configuration (see figure[T]i. 

Definition 1 F has the property Spot[x, y, I, r] for x £ Z and y,l,r £ N if 

• Fy c is a bijection for all configurations c; and 

• FJj c is a constant function for all c and all z £ [x — I; x + r] \ {x}. 

thus fulfills Spotfa;, y, I, r] if and only if, for any z £ [x—l;x + r],B z C B y is equivalent to x = z 
(considering that "B z C B y " is a false statement when B z is undefined). 

Lemma 1 IfF has the property Spot[a;, y, I, r], then F V (Q Z ) contains all the words of size max(L r) + 1. 

Proof: Let us suppose, without loss of generality, that I is no larger than r. Let q = (go, • • • , Qr) £ Q r+1 - 
We are going to construct c £ Q z such that i 7 ' !/ (c)o ! ....r = <?■ Start with an arbitrary c £ Q z . We can 
first modify c r _ x in such way that F y (c) r = q r ; then we can change c r _ x _i, on which F v (c) r does not 
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Fig. 1: Property Spot [a;, y, I, r\. Gray zones correspond to cells whose state does not change (either fixed in the 
initial configuration or kept constant in the y-th iteration of the CA). 

depend, so that F y {c) r -i = q r -i', and so on, until we choose c- x , on which F(c)i im- r does not depend, 
so that F(c) = q . □ 

For instance, (J) fulfills Spot[0, 1, +oo, 0], which implies that it must be surjective. 

Lemma 2 Let N v be the neighborhood of F v . If F fulfills Spot [a:, y, I, r] and Spot[x', y',l' , r'] with 
[x' - V) x' + r'] + M y ' C [-/; r], then Spot[ir + x' , y + y' , I', r'] also holds. 

Proof: By definition of the neighborhood, [x' — I'; x' + r'] +J\f y ' C [—l;r] implies that F y (c)u/_// ;a ./ +r /] 
is a function of Cr_; ;r i. Applying that to c = <J- X o F y (d), we get that the restriction of F y+y (d) to 
[x + x' — I'; x + x' + r'} is a function of F y (d)i x _i ;x+r \. Thus we get: 

• F y+y \d) [x+x ,_ v . x+x , +r , ] \ {x+x ,} does not depend on d ; 

• F y+V (d) x +x' depends only on F y (d)[ x _i ;x+r ], which in turn, according to Spot [a;, y,l,r], depends 
injectively on do. 

□ 

The central idea of the paper is to study the set of parameters (x,y,l,r) for which the property 
Spotfa;, y, I, r] holds, and use that set to obtain necessary conditions for simulations between cellular 
automata. However, we won't use the set of parameters directly because the simulation relation is invari- 
ant by space-time rescalings and this set is not. Instead we will look at 'scale-free' structures inside this set 
of parameters. More precisely, given some integer p, we look for infinite geometric progressions of order 
p in the set of parameters. Hence we obtain a kind of fingerprint for each CA which is well-behaved with 
respect to space-time transformations involved in the simulation relation (Theorem [T] below). Moreover, 
as shown by examples developed latter in this paper, this fingerprint is closely related to the self-similar 
structure observed in typical space-time of some linear CA. Technically, this is how the definition goes. 

Definition 2 For a CA F and an integer p > 2, we denote X p (F) the set of points (x, j)elx [0; +oo) 
such that for some k £ N,for every large enough n £ N, F fulfills Spot[a;p™, yp n , p n ~ k , p n ~ k }. 

^2(©) is for instance the set of points (x, j/)eRx [0; +oo) that can be written x = and y = 
with a, b e N and B a C _B b : its restriction to R x [0; 1] is the dyadic part of a (shifted) Sierpihski triangle. 

It can be noted that X p is necessarily of measure 0, and is self-similar, since by Lemma [2] every point 
of X p is the tip of a small copy of X p within itself. One can also notice that if F is not surjective, 
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then X p (F) is reduced to the singleton {(0, 0)}. Indeed, if X P (F) is not reduced to a singleton, then 
according to Lemma[TJ the image of F contains every finite word, which implies, by compactness, that F 
is surjective. 

We now detail, in a series of properties, how X p is modified under the action of the transformations 
involved in the simulation of a CA by another. First, the shift. Let s z be the transformation of the plane 
defined by S z (x, y) = (x + zy, y). The following property is obvious, by definition of X p . 

Property 1 X p (a z °F)=s z (X P (F)). 

Let us now consider iteration and grouping. Let g f be the transformation of the plane defined by 
Qt(x, y) = (x, |): notice that q p (X p (F)) — X p (F). Let f m be the transformation of the plane defined 

by f m (x,y) = (^,y). 

Property 2 X p (F t ) is a dense subset ofQ t (X p (F)). 

Proof: The inclusion is immediate from the definition. What might be slightly less immediate is why 
these sets are not obviously equal. Given the definition, X p (F) must be included in M. p , where M. p is the 
set of real numbers having finite p-adic expansion. Actually, we do have g t (X p (F)) D M. p = X p (F t ), so 
the equality without the intersection is certainly true if t £ M p , not quite so in general. Let us now prove 
the density. 

Let (x, y) £ X p (F). We want to find a sequence (x n , y n ) of points of X p (F) converging to (x, y) such 

that for all n, y n £ tM. p . For a finite sequence of integers = i W) o < Vi < ■ ■ ■ < hi.i (I is a constant 

l 

independent of n to be fixed later), we define r\ n = p~ l "- j and (x ni y n ) = i] n (x, y). We have three 

3=0 

requirements: 

• {%m Vn) must converge to [x, y): it is sufficient to have lim i n i = +oo 

n— >+oo ' 

• (x n ,yn) must be an element of X p (F). This is guaranteed as long as i n ,j+i — in,j is always 
large enough. More precisely, by definition of X p there exists some k such that for every large 
enough integer n, F fulfills Spot[xp n , yp n ,p n ~ k ,p n ~ k ]. Therefore, if i n j — i n j-i is large enough 
(depending on k and the neighborhood of F), we get from Lemma|2]that (1 + p~ ln > l+ ' ln < l -' L )[x, y) 
is in X p (F). By recursion on I, we get ultimately (x n ,y n ) £ X p (F). 

i 

• y n must be in tR p , which means the integer p 1 ™' 1 p~ ln -j must be a multiple of t. 

3=0 

So, it all boils down to finding increasing integer sequences = Iq < i\ < ■ ■ ■ < i\ where ij + i — ij is 

i 

arbitrary large, and such that t divides p %l P~ ij ■ That is clearly possible: the sequence of powers of p 

is ultimately periodic modulo t, so if we choose the ij-s spaced by multiples of this period and I = t, we 
can easily meet the conditions. □ 

Property 3 X p {b m o F o is a dense subset off m (X p (F)). 
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Proof: Let G = b m o F o Let x £ Z and y, Z,r e N with i > 1 and r > 1. First, it fol- 

lows from definitions that, for any configuration c of F, G v x b is constant if and only if, for all 
z G {mx — m + 1, • • • , mx + m — 1}, F^ c is constant. Moreover G v x b , c -. bijective implies F^ c bi- 
jective. This shows that if G has property Spot[x, y, I, r] then F has property Spot[ma;, y, ml, mr]. 

Now suppose that F has property Spot[mx, y, ml, mr] and fix some configuration c of F. Then it is 
straightforward to check that G v x b ,n is bijective (because it sends each component of Q rn on itself) and 
G v z b is constant for any z e [x — I; x + r] \ {x}. 

We have shown that F has property Spot[mx, y, ml, mr] if and only if G has property Spot[x, y, I, r]. 
Thus we have 

(x,y)eX p (G) ^ (mx,y) E X p (F). 

This implies X p (G) C f m (X p (F)). To prove the density, it is sufficient to prove that X p (F) n ml p is 
dense in X p (F) which can be done using the same argument as in the proof of property [2] □ 

It only remains to consider the case of the sub-automaton. 

Property 4 If G C F then X P (F) C X P (G). 

Proof: It is straightforward to check that if F has property Spot[x, y, I, r] then so does G. The property 
follows. □ 

Properties[T]|2]|3]and|4]prove the following theorem. 

Theorem 1 If F simulates G, then there exist rational numbers {3 and a, 7 > such that for every integer 
P>% K a ,0,-y(X p (F)) C X p (G), where ir a ,p,-f(x, y) = {ax + /3y, jy). 

The determination of X p is not easy in general, but the following basic facts can be established straight- 
forwardly from the definitions: 

• if F is a shift, then X p {F) is a line passing through the origin; 

• if i 7, is nilpotent (i.e. 3t s.t. F f is a constant function), then X p (F) — {(0, 0)}; 

• X P (F x G) = X P {F) n X P (G). 

Theorem [T] above shows that X p (F) represent obstructions for F to simulate other CA: the bigger 
X p (F) is, the smaller the family of CA F can simulate. Using the basic facts above, we can give some 
concrete formulations of this intuition. 

Corollary 1 Let p > 2 be an integer and F a CA. Then we have: 

• If F simulates the identity, then X p (F) must be included in a line passing through the origin; 

• If F is intrinsically universal, then X p (F) = {(0, 0)}; 

• If F is reversible universal (i.e. it can simulate any reversible CA), then X p (F) = {(0, 0)}; 

Proof: All items use Theorem [T] Item 1 and 2 are direct consequences of the computation of X p for the 
identity and nilpotent CA (an intrinsically universal CA must simulate any nilpotent CA). Item 3 uses the 
fact that a reversible universal CA must simulate a x cr -1 whose X p is a singleton. □ 

The purpose of the next section is to focus on a class of CA that generally have more interesting X p : 
linear cellular automata. 
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Fig. 2: Spacetime diagram of up to a large power of 2. Also X2{&). Time goes from bottom to top 

More often than not, one can get a good idea about what X p looks like just by examining the spacetime 
diagram. We think in particular of linear CA in the sense of IG NWlOl . In this case, Q = R d , where R is 
a finite abelian ring, and d some positive integer. The algebra of CA that are homomorphisms of (R d ) 
is then isomorphic to ^d(R) [u, u^ 1 ]: read section 1 of IG NWlOl for details. 

If F is such a linear CA and if denotes the neutral element of R d , the sets X p can be derived from 
the functions F^- where denotes the uniform configuration everywhere equal to 0. Indeed, for any 
configuration c, we have: 

F% c bijective (resp. constant) -4=4> F^ bijective (resp. constant) 

In the sequel we denote F^- by F%. The remainder of this section focuses on reversible cellular 
automata. 

3.1 0: a reversible CA which cannot simulate the identity 

Let us look at a more interesting example. The alphabet is now (Z2) 2 , and the transition is given by 

1 



6 = 



1 + l 



Since it already serves as a red thread through [GNW1Q], we will pass very quickly on it. Let us notice 

r 1 + 1 + u l\ 



here that, since its determinant is 1, it is reversible, and that its inverse is O 1 = 



1 
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Obviously, 9 simulates its own inverse: in fact 9 1 C v 9 with ip = f ^ J 

Figure [2] represents the spacetime diagram of 6 up to a large power of 2, for an initial configuration 
consisting of one single nonzero cell. 9 is "well-behaved" in the sense that these spacetime diagrams, for 
increasingly large powers of 2, converge to X2(Q)- It thus gives in a sense a purely visual proof of the fact 
that 9 does not simulate the identity. Of course, this requires actually some background knowledge, in 
order for the proof to be correct. One must know that 9 is a linear CA, and that X p actually corresponds 
to its limit spacetime diagram, or at least is not limited to one line. While X p is not defined in IGNW10I . 
the information given there on the way to describe the limit spacetime diagram by means of a substitution 
system justifies this assertion. The crucial point is that any block that is not empty contains a reduced copy 
of the whole pattern, which means that in the neighborhood of any non-white point in the limit spacetime 
diagram, there is a copy of the whole thing, whose tip is then a point in ^(9); therefore ^(9) is dense 
in this pattern. And so, adding that 9 is its own mirror image, we get: 

Proposition 1 9 simulates its mirror, its inverse and its dual, but cannot simulate the identity. 

3.2 r : a life in pictures 

Let us now provide an example of a CA that is both space- and time-asymmetric, in the sense that it cannot 
simulate any of the CA derived from it by inverting space and/or time. This will be 

1 

1 u | e ^3(Z 2 )[w,m _1 ] 

1 u 



Its inverse is given by T 1 



We are going to give only the proof that it does not 



ii 2 ii 1 
u 1 
1 

simulate its inverse: the proof of the two other results would add only length to this article, and can surely 
be left as an exercise to the reader. 




(a) X 2 (F) (b) XziF- 1 ) 

Fig. 3: X2 with the second coordinate restricted to [0, 1] (time goes from bottom to top) 
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Let us imagine for one blissful moment that we know X 2 to be accurately represented by Figures 3(a) 



and 3(b) (actually these figures are mirror images of spacetime diagrams up to a large power of 2). How 
do we conclude then? 

Supposing that V simulates T _1 , we know from theorem[T|that for some a, j3, 7, ir a ,p,-y (^(r)) should 
be included in ^(r -1 ). Since there are only two lines passing through the origin in ^(r -1 ), 7r Q ( 3 7 
must send respectively E(0, 1) and R(— 1, 1) on M(0, 1) and K(-2, 1), which implies (3 = 0. Now if we 
consider the lines joining these two axes, they have slope | for T, 1 in r _1 , which means a = 27. So, 
if r simulates its inverse, X 2 (T) should be, modulo a change of scale, included into X 2 (r~ 1 ), which is 
clearly not the case. 

To make this proof rigorous, we need a tool to prove properties of X 2 for T and r -1 . We are going to 
follow section 3 of [GNWIO], which gives a procedure to derive, from the transition matrix of the CA, a 
substitution system generating the Green functions (see Proposition 4 of |GNW10|). More precisely, we 
will associate to each CA F a 2 x 2 substitution system, that is a finite set E and a function e : Z x N ^ E 
such that: 

• F y is a function of e(x, y)\ 

• for i, j £ {0, 1}, e(2x + i, 2y + j) is a function of e(x, y) and i and j. 



The next two subsections give the substitution systems for T and T , and subsection 3.2.3 uses them 
to formally prove negative result concerning simulation. 

3.2.1 A substitution system for T 

The minimal polynomial of T is X 3 + X 2 + (1 + u 2 )X + 1, so we have the following recurrence relation. 

Vx £ ZVn, y £ N y<3-2 n =^ T^+y = Tf +1+y + Vf +y + + Y v x (1) 

Now we define Cij(x, y) in the following way: these are the coefficients in Z/2Z such that for every 
function (x, y) H> S| fulfilling equation ([T]l in lieu of T, 

2 

For every x £ Z, y £ N and s,t £ {0, 1}, we have 

ZllXl =Y,a (x-i, y)El l+s + ai (x- i, y)E 2 2 +l s + a 2 (x - i, y^t+s- 0) 
In the case s = t = 0, we have the following derivation: 



-2x 



J2 ot Q {x - i, y)z? 2l + at(x - i, y)H| j + a 2 { 



2i 



J2 a Q {x - i, y)H^ +a t (x- i, y)E% + a 2 (x - i, y) {z? 2l _ 2 + 3^_ 2 + S^) 



J2 (a (x -i,y) + a 2 (x - i, y)) 3^ + a 2 (x - 1 - i, y)3. 



(4) 



+ (ai(x -i,y) + a2(x - 1 - i,y))E 



■='2 

2i 
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which is to be compared with the definition of aj : 



i j=0 



(5) 



The comparison shows that aj(2x, 2y) is a function of ctj(x — i, y) for some values of i. V is peculiar 
in that aj(2x + 1, 2y) = 0, which simplifies our work. The same operation now has to be performed for 

otj(2x,2y + l). 

s 2x +1 = T, a o{x-i,y)'Bl i + a 1 El i + a2(x-i,y)El i 

i 

= Y< a o{x-i, y)^\i +ai{x- i, y) (S^ + + 5^_ 2 + S^) 



+a- 2 (x - i, y) (S^ + S^_ 2 + + -21-2) 

{ai{x -i,y) + a 2 (x - 1 - i, y)) 

+(a- (x - i, j/) + ai(a; - 1 - i,y) + ai(x -i,y) + a 2 (x - 2-i,y) 
+a 2 (x -l-i,y)+a 2 {x-i, 2/))E 2i + a x (x - i, y)~? 2l 



Using the representation 



a 2 
ax 
a 



we get the following substitution. 



a.(x,y) 



a.(2x,2y + l) 


a.(2x + l,2y + l) 


a.(2x,2y) 


a.(2x + l,2y) 


11 



Qi(x,y) 







a (x, y) + ui (x - 1, y) + «i (x, y) + a 2 (x - 


2, 2/) + "2 (z - 1, y) + a 2 (x, y) 





ai{x,y) + a 2 (x - 







ai(x,y) + a 2 (x - 


i,y) 





oti{x - l,y) 







ao (a;,y) + a 2 (a; 


,y) 






(6) 



This needs some grouping; for instance, in the present situation, the substitution scheme uses ot\{x — 
1,2/), which is not an information contained in the initial cell. For instance, if we want to determine 
ao(2x,2y + t) for t € {0,1}, we need to know a (x 7 y), a\(x,y), a 2 (x,y) and 0:2(2; — 1,2/)- Th e 
smallest grouping that will allow us to carry all that information is 



a 2 (x-3,y) a 2 {x- 


2,2/) 


a 2 (x - 


1,2/) 


a 2 (x,y) 


a\(x — 


2,2/) 


oei(x — 


1,2/) 


ai(x,y) 






a (x - 


1,2/) 


ao(x,y) 



This gives us an alphabet of size 2 9 = 512, and the substitution scheme is 
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a b c 


d 


e f 


9 


h 


i 



/ g 

a+b+c+e+f+h b+c+d+f+g+i 
c + g 


/ «? 

b+c+d+f+g+i 

c + .g 


b+f c+g 
b c 
d + i 


b+f c+g 
c 
d + i 





abed 

e / g 
h i 


The initial state for this substitution system is 





D 





■ ■■ , and to a cell 
















d + 


i 


c 




position (x, y) corresponds the Green function T% = 


■( 


c 




b + d + i + q c 


l > \ 






9 




c + f b + 


d + i J 



For a letter x in {a, 6, . . . , i}, let X denote the cell where x has the value 1 whereas all other letters 



are set to 0. We can notice that A, E and H are completely equivalent: they all substitute to 



E 












and project onto in the computation of r*. We can therefore simplify this system a bit by putting 

A = E = H = 0: 

Whereas we have a theoretical number of 2 5 = 32 different states in the substitution scheme, only 1 1 
of them are accessible from the initial state, namely plus the ones represented in Figure [5] This graph 
has two strongly connected components, one composed of BD alone, the other of the remaining vertices. 
In particular, from any state of the substitution system that has been accessed from the initial state and 
that is neither nor BD, there is a path to D; therefore there must be a point of X2 in the corresponding 
square. 

3.2.2 A substitution system for r 1 

We now have to perform the equivalent analysis for which we will name f2, in order to avoid possible 
confusions with negative exponents. The minimal polynomial of £1 is X 3 + (1 + u 2 )X 2 + X + 1, so now 
the recurrence relation is 



Vx e ZVn, y E N y < 3 • 2" => n^+v = fif +1+y + nf^Vi + ^T +v + (7) 
We introduce j3, which is to il what a was to T in Section [3.2.1| 



= £ A>(z - i, y)S' 2l+s + Pi(x - i, y)S™ + /3 2 (x - i, y)E i +l 



(8) 
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Fig. 4: Fifth step of T's substitution system (time goes from bottom to top). 
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Fig. 5: Transition graph of the substitution system: an arrow from state si to state S2 means that S2 can be obtained 
after a finite number of iterations starting from s\ . 

We then get the following decompositions. 

s 2z = J2/ 3 o{x~i,y)^ 2l + P 1 {x-i,y)El l +f3 2 (x-i,y)Ej i 

= jzPo{x-i, y)S° 2i + /3 1 (x- i, y)~ 2 2l + fc{x - i, y) (E 2 2i _ 4 + ~ l 2l _ 2 + E%_ 2 + 3^) 

i 

= E (Po(x-i,y) +fc(x- 1 - i,y)+p 2 (x - i, y)) + f3 2 (x - 1 - i,y)E\ t 
+ (Pi(x -i,y)+0 2 (x-2-i, y)) S| 

This is to be compared to this definition of j3j : 

2 



— 



Likewise, for y i-> 2y + 1, we get 



= E (Pi (x-i,y)+h(x- 2-i,y))Z% 



The minimal grouping is now. . 



(9) 



£^&(2:r-i,2y)S| (10) 

i j=0 



+\fio(x-i,y)+Pi(x-i,y) 

+/3 2 (x -2-i,y) + f3 2 (x - 1 - i,y) + (3 2 (x - i, y)\E} 2l ( n ) 
+[/3 x (x- l-i,y) + px(x-i,y) 

+/3 2 (x - 3 - *, y) + fo{x -2-i,y) + (3 2 {x -1-i, y)]3§ 4 
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/3 2 (x-5,y) f3 2 (x-4,y) fa(x 


-3,2/) 




-2,2/) 


P 2 (x- 




Pi{x,y) 


(3i(x 


-3,») 


Pi(x 


-2,2/) 


(3i(x- 


1,2/) 


Pi{x,y) 










A)(+ - 


1,2/) 


Pa{x,y) 



. and the corresponding substitution scheme is given by 



a+b+c+g+h b + c + d+h + i c+d+c+i+j 
c + d + c + , + k d + e + / + j + I 
d + j 


a + b + c + 9 + h b + c + d+h + 1 c + d + e + i + j 
c+d + e + t + k d + e + / + j + ! 
d + j 


(] b + h (I c + i d + j 
d e 

e + / + 1 


b + b c + i d + 3 11 
d e (1 
e + / + I 11 



The initial state is ••• L K ••• , and (T 1 ) y is given by 



I + h + f + d + b i + e + c j + d 
i+e+c l+j+f+d e 
j+d e l+f 

A being equivalent to G, B to H and F to L, we get the simpler 



c d 


e 


g h 


i j 




k I 



c+d+h+i 
d + e + g + h + i + k d + e + j + I 

c+e+i 


c+d+h+i c+d+e+i+j 
d+e+g+h+i+k d + e + j + Z 

d + j 


c + i 

d + ft e 

d+e+j+l 


c+i d+j 
d+ft e 
e + l 



This first simplification makes G equivalent to K, C to IK and J to Z?i7, so we finally get 



d 


e 


ft 


i 


fc 


Z 



e+h+i+l 
h+i+k+l 
e + i 


d+e+i 
e+h+l 
d+k 


e + i 
d+e+h 
d+e+l 


d 
e + i 
d+e+ft+i+Z 



. which results after five steps in Figure [6] 
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Fig. 6: Fifth step of F 



1 's substitution system (time goes from left to right). 
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(0, 1) 



(0,2) 




4 2 

3 ' 3 



*2<r) 



-l 



Fig. 7: Partial knowledge about ^(r) and ^(T -1 ). Points in black are known to belong to the set while points 
in gray are known not to belong to the set. Remember that X 2 is invariant by homothetic transformations of center 
(0, 0) and factor 2 l (with i any integer). 

3. 2. 3 Final arguments 

It now remains to be proven that Figure[3]does represent X 2 for T and As such, this does not mean 
much; actually, we need to prove a few features of X 2 that would suffice in order to conclude that T does 
not simulate Namely, we want to justify this series of assertions: 

(i) X 2 {T) contains the (half-)lines R+(0, 1) and E+(l, 1). 

(ii) X 2 {T) contains the segment [(0, 1); (§, §)]. 

(hi) No point of X 2 (T~ 1 ) lies in the interior of the triangle with vertices (0, 1), (|, |) and (|, |). 

This is enough to conclude, because (iii) implies that the only possible half-lines starting at the origin 
and included in ^(T -1 ) are the vertical axis and that of slope ^; and the segments joining these lines, 
if they exist, must have slope —1. Therefore it is impossible to send X 2 (T) into X 2 (T~ 1 ) by a 7r a ,p,-y 
transformation (see figure |7]i and Theorem [l] concludes. Since X 2 (T) is just the symmetric of X 2 (T~ 1 ) 
with respect to the vertical axis passing through (0,0), the same reasoning with Theorem [T] shows that T 
cannot simulate T. 

Proposition 2 T can neither ^-simulate its inverse Y^ 1 nor its dual T. 

We now prove the three assertions above successively using the substitution systems derived earlier. 
Property 5 X 2 (T) contains the (half-)lines M+(0, 1) and 1). 

Proof: First, by looking at the images of D, G, B and F by the substitution system of T, we prove by 
recurrence that: 

• T(0, n) = D if n is even and G else, and 

• T(n, n) = B if n > 1 is even and F else, 

where T is the fixed-point of the substitution. We deduce from the former observation that every letter in 
the substitution system, except for BD, contains a point in X 2 , that the (half-)lines M+(0, 1) and R+(l, 1) 
are in X 2 (T). □ 
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Property 6 X 2 (T) contains the segment [(0, 1); (|, |)]. 
Proof: The substitution system of T is such that: 



BD 



BDG 









BD 






G 


CF 


B 


C 



CF 
G 



BD 





BDG 


CF 



G 


CF 


D 


C 



We deduce that any pattern of the form 



where X is either BDG or G, is sent to a pattern of the form 



Now, observing Figure (jij, one can see a discrete segment of slope — | made of the above pattern starting 
from the top-left position and reaching the upper-diagonal. Since we know that all the cells appearing on 
this discrete segment, namely G, CF and BDG (plus an end point that is, depending on the parity of the 
scale, B or F), contain a point of X2, it just remains to show by recurrence that a segment of this form is 
present at every scale, which is immediate. □ 



Property 7 No point 0fX2.iT 1 ) lies in the interior of the triangle of vertices (0, h), (i, g) and (4 



^3> 3' 



Proof: By induction, we can prove that depending on the parity of the step, this triangle takes alternatively 
the forms presented in Figures [8] and [9] which represent the corresponding triangle in the substitution 
system, supposing the pair of initial blocks represents a rectangle of height 1. For instance, Figure [6] 
showing the fifth step, exhibits in this position a triangle of the form presented in Figure [8] 

The proof that each of the figures substitutes into the other one is purely mechanical, and essentially 
done by the very existence of Figure [6] where the first five steps of substitution are readable. □ 



4 Discussion 



We gave a new necessary condition for the simulation of CA and applied it to solve a few open questions 
of the form 'Does there exist a reversible CA that simulates such and such but not such and such?'. 
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DHL 
DHL 



DHL 
DHL 
DHL 
DHL 
DHL 



D 



Fig. 8: Odd steps (time goes from bottom to top). 



Noticeably, each time we were able to answer this question, it was in the positive, which is one general 
reason why we would expect the same answer for other closely related questions of the same sort that 
remain open. 

Our method is tailored to be applied to linear CA. Their practical advantage is that much of the infor- 
mation is present in their spacetime diagram, and therefore easy to access and comprehend. For instance, 
with our theorem in mind, a blink at Figure [2] is enough to suspect that O cannot simulate the identity. 
It then remains to check rigorously that the pattern does represent X^ accurately, but that part is purely 
mechanical, if a bit tedious. Let us now finish with two questions. 

Why did the authors resort to a 3 x 3 matrix? Couldn't they find anything simpler? No, they could not. 
Actually they conjecture that every 2x2 matrix simulates its inverse, which interestingly enough reduces 
to deciding whether every matrix simulates its transpose. 

Does there exist a CA that can simulate the identity, but not its inverse/dual? The correct answer is 
'probably, and V x id is a good candidate' . However, our theorem is not really helpful in this case, since 
the Xp-s of this CA are trivial. Hopefully some hybrid can be created by merging it with [DMOTllb 
theorem 3.4] and made available to the masses in the future. 
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Fig. 9: Even steps (time goes from bottom to top). 
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